Impurity modes and effect of clustering in diluted semiconductor alloys 
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The variation of TO zone-center vibration spectra with concentration in mixed zincblende-type 
semiconductors can be understood within a paradigm of unified "one bond - two modes" ap- 
proach, which has been recently outlined as a rather general concept,^ and emerges from a 
number of previous experimental and theoretical studies. ^'^ The crucial issue is that the vibra- 
tion frequency, associated with a certain cation-anion bond, depends on the length of the latter, 
and the bond length, in its turn, depends not only on the average alloy concentration, but on 
local variations of it. In an (A,B)C substitutional alloy, the A-C bond length differ in A-rich 
and A-poor regions, yielding a splitting of the A-C vibration frequency. Such splittings can be 
measured and reproduced in first-principles calculations. 

An analysis of vibration spectra helps to get an insight into the structural short-range (clus- 
tering) and long-range (formation of extended chains of certain cation-anion pairs and other 
structural motives at the mesoscopic scale) tendencies. For this however, one needs first- 
principles benchmark calculations for representative model systems (see, e.g., Ref. 4 for the 
ZnSe-BeSe alloys). The simplest yet important result from first-principles calculations is a pre- 
diction of how the impurity phonon mode evolves as isolated (distant) impurities get clustered. 

In the present contribution, we outline the results of first-principles calculations of phonon 
frequencies and vibration patterns, in the dilution limits of several mixed semiconductor al- 
loys, Bea;Zni_a;Se, Gaa;Ini_ajAs and Gaa^Ini-^^P. The calculations have been done by the Siesta 
method^ for cubic 64-atom supercells, with one or two cation atoms substituted by impurity 
species, that corresponds to impurity contents of 3% and 6%, respectively. The initial uncon- 
strained structure relaxation for each supercell chosen was followed by a calculation of phonons 
by finite displacement technique. Each of the atoms in the supercell were subject to 6 consecu- 
tive small cartesian displacements, and the forces induced on all atoms in the supercell resulted 
in corresponding force constants. 

The analysis of results is the simplest for Bea;Zni_a;Se, a system with large contrast in 
masses and elastic properties between its parent compounds, that leads to a big separation be- 
tween Zn-related and Be-related phonon modes (see Ref. 4 for details). Fig. 1 shows the phonon 




Fig. 1. Phonon density of states 
of Be, after first-principles cal- 
culations for BenZn32-nSe32 
supercells, with 1,..4 Be 
cations neighbouring the same 
Se anion. The supercell for 
n=4 is shown in the inset. 
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density of states for Be, for the cases of single impurity and impurity pair (both Be atoms being 
neighbours to the same Se). Moreover, the cases of three and four Be atoms grouping around 
the same Se are considered for this system (the Be4Zn28Se32 supercell is shown in the inset to 
Fig. 1). One clearly sees that the single vibration peak of an isolated impurity (a nearly triply 
degenerate mode of Be in tetrahedral Se4 cage) splits into two peaks, separated by ^40 cm~^, 
as another Be impurity is introduced. The origin for this splitting is a diversification of Be-Se 
bond lengths, as two Be share the same Se atom. Namely, the bonds to outer Se atoms of the 
Be-Se-Be cluster can be much more efficiently shortened in the course of structure relaxation 
than the bonds to the central Se. Shorter bonds imply larger force constants and hence higher 
vibration frequency (again, see Ref. 4 for a more detailed discussion). It is noteworthy that an 
addition of further Be atoms as neighbors to the central Se brings in some additional structure, 
but leaves in place the initial major splitting into two big groups of peaks, those characterizing 
a "Be-poor" environment (at ~480 cm~^) and "Be-rich" environment (at ~440 cm"^). 

In the following we compare this initial splitting of single-impurity mode into a doublet 
of interacting-impurities modes with the experimentally observed onset of the "l-bond^2- 
mode" behaviour, discussed previously. Notably, the impurity mode at ^440 cm~^ in the 
x^O limit of Bea,Zni_a;Se develops into two branches, at ^^450 cm"^ and ~410 cm"^, at small 
values of X - see Fig. 1 of Ref. 6. It should be noted (and was already discussed in Ref. 4) 
that the absolute values frequencies are shifted upwards in our calculation, with respect to their 
experiment values, due to a slight overbinding caused by the local density approximation to 
the exchange-correlation. However, the frequency difference between the single-impurity and 




Fig. 2. Simplified l-bond^2- 
mode TO (thick lines) percola- 
tion schemes of (InGa)As and 
(InGa)P. The generic fraction 
of bonds corresponding to each 
TO branch is indicated in the 
bottom panel. The bond-related 
optical bands simply obtained 
by linear convergence of the 
parent TO and LO frequencies 
onto the related impurity fre- 
quencies are shown as shaded 
areas. 



double-impurity modes is fairly reproduced. 

The general behaviour of phonon modes depending on the Be concentration, "streamlining" 
the experimental details from Fig. 1 of Ref. 6, is shown and discussed in Ref. 7. It was argued^ 
that the similar schema holds for other semiconductor alloys, even if the position of branches 
and their splittings may differ. Fig. 2 reproduces a part of experimental findings schematically 
presented in Fig. 4 of Ref. 1. Below we discuss the ab initio predictions of the single impurity 
vs. impurity pair -related phonon behaviour for (Ga,In)As and (Ga,In)P alloys. 

For the identification of zone-center phonon modes it is convenient to refer to q-projected 
phonon density of states (PhDOS), which can be introduced as 

hiu, q) = E| E exp(qR„) ~ 

where Af{(jj) is phonon eigenvector for atom a at frequency cj, and - an arbitrarily chosen 
group of atoms (say, those of a given chemical species). The q=0 -projected PhDOS of a 
host atom selects the zone-center vibrations which contribute to the Raman spectra, like e.g. 
a TO-mode of bulk InAs at ^--^220 cm~^ is prominent in the PhDOS of In in Fig. 3. For Ga, 
which is a single impurity, the q-projected PhDOS coincides with the full vibrational density 
of states and spans the whole interval of frequencies, from to 250 cm"^. Yet, not all of these 
modes are zone-center-like, if one takes into account the vibration of the As sublattice. A direct 
inspection of different vibration patterns reveals, as genuine zone-center vibration modes, a 
triplet at 246 cm"^. In the q=0 -projected PhDOS of In (Fig. 3) and of As (not shown, as it is 
practically identical with that of In), a split-off peak on the high-frequency side of the main TO 
line does also reveal the presence of this impurity mode, which is an in-phase vibration of Ga 
impurity with the In sublattice, in counter-phase to the As sublattice. 

Turning to the analysis of impurity pair, we find several modes, which correspond to in- 
phase vibration of Ga atoms with the In sublattice and in counter-phase with the As sublattice. 
These modes span the range of 248-251 cm"^. 

On the opposite end of the concentration scale, the vibration of isolated In in GaAs, as a 
heavy impurity, mixes up with a "continuum" of bulk modes. Still, the major impurity modes 
can be identified near 238 cm"^. The impurity pairs vibrate in phase between themselves and 
with the cation sublattice in the region 235-241 cm"^, again in agreement with the measured 
vibration zone-center frequencies shown in 
Fig. 2. The first-principles calculations yield 
therefore, on both low-concentration sides 
of the concentration scale of (Ga,In)As, 
the zero splitting from single-impurity 
frequency to the two-mode regime. This is 
consistent with experimental observations 
summarized in Fig. 2. 



Fig. 3. q=0 -projected phonon density of states 
of Ga and In in Gailnai AS32 supercell. 
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Fig. 4. Two in-phase 
vibrations of Ga impurity 
pair with strong zone- 
center contribution in the 
vibrations of the anion 
sublattice, according 
to calculations in the 
Ga2ln3oP32 supercell. 

For (Ga,In)P the situation is different. In the experiment, the two-mode regime which sets 
on under the Ga doping is characterized by a quite large sphtting (of ^20 cm"^) between the 
two Ga-P modes. Indeed we could recover this splitting in the calculation. Fig. 4 shows the 
vibration patterns of two maximally split in-phase double-impurity modes. The softest of these 
modes (322 cm"^) is a "longitudinal" displacement of two Ga atoms against the P atom in be- 
tween, making one Ga-P bond longer at the expense of the other. The hardest mode (352 cm"^) 
involves a different "two Ga against P" displacement, which lies in the Ga-P-Ga plane as in the 
first case, but occurs in the perpendicular direction, bending the Ga-P-Ga angle. This vibration 
is slightly harder than the single-Ga impurity mode whose calculated frequency is 350 cm"^. 
Taken together, the arrangement of calculated vibration modes agrees well with the experi- 
mental observations summarized in Fig. 2 (up to a small systematic blue shift of calculated 
frequencies with respect to experiment, discussed above). 

On the Ga-rich side of the schema in the lower panel of Fig. 2, no noticeable splitting was 
detected in experiment, and again this is consistent with the calculation: single-In impurity 
modes at 360 cm"^ develop, on addition of a second In impurity, into a bunch of two-impurity 
modes, which are confined in a quite narrow range, 357 to 362 cm"^. 

Summarizing, we have shown, on discussing impurity limits in three semiconductor alloys, 
how the important parameters of the concentration-frequency diagram for mixed semiconduc- 
tors discussed in Ref. 1, the single impurity frequencies and their initial splitting on the onset 
of one-bond^two-mode behaviour can be extracted from first-principles calculations. 

We acknowledge the access to the calculation resources at the CINES in Montpellier (projet 
N° pli2623). 
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